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We consider the quantum correlations, i.e. the entanglement, between two systems uniformly 
accelerated with identical acceleration a in opposite Rindler quadrants which have reached thermal 
equilibrium with the Unruh heat bath. To this end we study an exactly soluble model consisting of 
two oscillators coupled to a massless scalar field in 1+1 dimensions. We find that for some values of 
the parameters the oscillators get entangled shortly after the moment of closest approach. Because 
of boost invariance there are an infinite set of pairs of positions where the oscillators are entangled. 
The maximal entanglement between the oscillators is found to be approximately 1.4 entanglement 
bits. 



I. INTRODUCTION 



The ground state of a quantum field is an extremely structured state. This is true even if the field is non interacting. 
A first indication that this is the case is that the propagator between two spatially separated points never vanishes, no 
matter how far apart the points are. That is correlations in vacuum extend over an infinite range. More subtle is that 
these correlations are such that in vacuum two spatially separated regions are entangled. This was first exhibited in 
T-H , Q by decomposing Minkowski space into two Rindler wedges, and quantizing the field in each wedge. It is also very 
\^ • closely related to phenomena like black hole evaporation or particle creation induced by cosmological expansion0,|^. 

It was later shown that the ground state exhibits quantum non locality: it is in principle possible to violate Bell 
^^O inequalities in the vacuum, both for free fields|3, 0, and for interacting fields'r'l, and that this property applies to 



> 



almost any quantum statep,Q- Recently B. Reznik considered a specific model in which two localized detectors are 
, coupled to the vacuum in space like separated regions and showed that the detectors can get entangled |l3 • In it 
• was shown that the correlations between the detectors could be used to violate a Bell inequality. These studies were 
CIh! later extended to the case of more than two localized detectors . 

1^ • In the present work we pursue the study of how uniformly accelerated detectors in opposite Rindler wedges get 
entangled. Contrary to previous work jlOj we shall study the case where the detectors are in equilibrium with the 
Unruh heat bath. Nevertheless, although both detectors have thermalised, we shall show that they get entangled for 
certain choices of parameters and for certain relative positions. 
^ , The model we shall study consists of two oscillators uniformly accelerated with identical acceleration a in opposite 
Rindler quadrants coupled to a massless field in 1+1 dimensional Minkowski space time, see Fig. ^for a depiction 



of the trajectories. This model is exactly soluble, which allows us to study the case where the oscillators are in 
equilibrium with the Unruh heat bath, and also to study the regime where the interaction between the oscillator 
and field is strong, and perturbation theory is no longer valid. We find (numerically) that, for certain values of the 
parameters, the two oscillators, both in thermal equilibrium with the Unruh heat bath, indeed get entangled. An 
interesting aspect concerns the position at which the oscillators are maximally entangled. Naively one would expect 
that this should occur when the squared invariant distance As^ = Ax^ — At^ between the oscillators is minimal (i.e. 
when the always spacelike separation between the two oscillators is minimal). But in fact, due to the dynamics of 
the oscillators which absorb and emit quanta at a characteristic rate, the entanglement occurs only for a short period 
slightly after the moment of closest approach. Note also that because of the boost invariance of the problem, there is 
an infinite set of pairs of points where the entanglement will be maximal. 

The model of an oscillator coupled to a massless field in 1+1 dimensions has already been extensively studied, 
both in the context of the Unruh effect [H El El El [13, 

and with the aim of understanding decoherence and 
thermahsation[Tl[l9ll20|. Much of our analysis is based on this earlier work. However during the present investigation 
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FIG. 1: Minkowski and Rindler coordinates. The Minkowski time and space coordinates are t and x; V,U — t it x are the 
Minkowski hght hke coordinates; tr, tl and pR, pL denote the Rindler time and space coordinates in the right and left 
Rindler quadrants respectively. Note that tr increases toward the future, whereas tl increases toward the past. The uniformly 
accelerated oscillators follow the trajectories pR = a~^, and pL = which are indicated in the figure by bold curves. Due to 
Lorentz invariance, there is an infinite set of pairs of points for which the squared invariant distance between the oscillators is 
minimal, such as the pair {Pl,Pr) or the pair {P'l,Pr) indicated on the figure. The two oscillators only are entangled shortly 
after having reached positions where As^ is minimal. For instance if the left oscillator is at Pl (or P'l), the right oscillator 
should be along its worldline, slightly after Pr (or P^,). 

we encountered a problem which had apparently not been noticed before, namely that the momentum of the uniformly 
accelerated oscillator has infinite fluctuations due to an infrared divergence. We show how these divergences can be 
controlled. 

Note that because the system we consider consists of two oscillators, we must use the tools which have developed for 
the study of entanglement of continuous variable systems, first considered in the seminal paper of Einstein, Podolsky 
and Rosen ■ The relevant tools will be reviewed below. 

II. THE MODEL 

We parametrize 1+1 dimensional Minkowski space as 

i = Pfl sinh a Tfl , x = pjiCOshaTR (\t\ < x) , 
t — ^pLSmhaTL , X = —pLCOshaTL {x < ~\t\) . (1) 

Let us note that th and — measure the proper time along the trajectories {pn — Pl — ct^^) of the two oscilators. 
Later on we shall make use of the null coordinates 

{\t\ < x) 
{x<-\t\) 

m < x) 

{x<-\t\) 

(2) 

Note that with these definitions tji, vj^ and ul increase toward the future, whereas r^, and un increase toward 
the past. 



V t + x = 

U = t-x= 



a expaufl 

^-1 , 



exp av L 



—a ^expaufl; 
a^^ exp auh 
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The problem of a single uniformly accelerated oscillator coupled to a massless field in 1 + 1 dimensions has been 
extensively studied Iisl ITgL IT^ . Here we generalise it to the case of two oscillators uniformly accelerated with 
the same acceleration a in opposite Rindler quadrants. The action describing this system is 

1, 



j dpRdTR5{pR- a^^) 
J dpLdTLS{pL - a^^) 



m f dqR 
dTR 

m f dqL 
dTL 



I dTR 



2 2 dqL . _i. 



(3) 



where qr and ql are the internal coordinates of the oscillators; and m, w, e are their mass, oscillation frequency, and 
coupling to the field. The momentum conjugate to the oscillator coordinates are 



PR = m- \-e(j){TR,a ) 

dTR^ 

PL = -m- \-e(p{TL,a ) 

dTL 



(4) 



The equations of motion are 



d^ qr 2 2 



'drl 



dqR ^ 
dTR 
_ d(f>{TR,a'^) 

dTR 
d(f>{TL,a''^) 



e^^S{pR - a ^) + e^^5{pL - a ^) 
dTL 



dTL 



(5) 



The first of these equations may be integrated to yield 



{X,t) = (t)°{x,t) + -qR{TRret) + T;lL{TLret) 



(6) 



where TR^Lrej, iare the values of tr^l at the intersection of the past light cone from {x, t) with the right and left 
trajectoriesj25j and where is the free field operator, solution of d^cfP — 9^0" = 0. We can reexpress the free 

(7) 



field solution as 



<tP{x,t) = (tP{t + x)+(tP{t~x) 



Inserting the solution for (p into the equation for qR^L yields the equations 



d^qR ^ dqR ^ 2 _ d(l)"{TR,a ^) 



m—^ + —-7— + muj qR 
dT^ 2 dTR 

d'^qL dqL 



m- 



drl 



2 dTL 

The equation for qR can be integrated to yield 



+ moj qL 



dTR 
i'°(TL,a-i) 



dTL 



(8) 



quiTR) = qo{TR) 



dTRxirR - T^je 



(9) 



A similar solution obtains for qL- Here x is the retarded propagator of and qo is solution ofmqo+e^qo/2+muj^qo — 0. 
The free solution go is exponentially damped. Henceforth we will neglect it. This corresponds to supposie that the 
oscillators were set into acceleration sufficiently far in the past and have reached thermal equilibrium with the Unruh 
heat bath. 
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Taking Fourier transforms, we can reexpress the solutions for the oscillators as 

TTl J 

m J 

-1 

X\ 



0°l(A) ^ I ^e-^-<A°(r^,a-i) . (10) 

Since the equations of motion are linear, these expressions represent both the solutions to the classical and Hcisen- 
berg equations of motion. In what follows we will take the initial state of the field cf) to be the vacuum state. The 
initial state of the oscillator is then irrelevant, since as we noted above, its state depends only on the state of the field. 

A key point in our analysis is that the oscillators are then in a Gaussian state. Indeed the initial state of the field, 
the vacuum, is a Gaussian state. And the internal coordinates qL,R and conjugate momenta pl.r of the oscillators 
depend linearly on the field operator. Therefore the oscillators are also in a Gaussian state. This means that the 
reduced density matrix of the oscillators, obtained by tracing over the field degrees of freedom, is entirely characterised 
by the expectation values of the first and second moments of the oscillator variables. In particular, as we will review 
in the next section, these moments completely characterise the entanglement between the two oscillators. 

It is immediate to obtain that the canonical variables have vanishing expectation value 

{qr) {ql) = {pr) = {pl) = . 
It is a more complicated task to compute the covariance matrix: 

({qR(TR)qR{TR)) {{qR{TR),pji{TR)})/2 {qR{TR)qL{TL)) {qR{TR)pL{TL)) \ 

{{qR{TR),PR{TB)})/2 {pR{TR)pn{TR)) {pRiTR)qL{TL)) {Pr{tr)Pl{tl)) , . 

{qdrMTR)) {qL{TL)pR{TR)) {qLirMrL)) {{qL{rL),PL{TL)}) /2 ^^^^ 

{pL{TL)qR{TR)) {pLiTL)pR{TR,)) {{qLiTL),PL{TL)})/2 {pl{tl)pl{tl)) I 

where {•, •} is the anticommutator. The covariance matrix depends on the positions tr , of the two oscillators. 
The expectation values along a single trajectory, such as {qR{TR,)qR{TR)) , are independent of the position along the 
trajectory, since we have supposed that the oscillators have reached a stationary state. The expectation values between 
operators on opposite trajectories, such as {qR{TR)qL{TL)) depend only on T — tr — tl, since by boost invariance it 
depends only on the invariant distance between the two oscillators : As^ = 4 cosha T/2. Thus ct is a function 
only of T. The detailed calculation of a will be carried out in section Hvl and in the appendix. 



III. ENTANGLEMENT IN CONTINUOUS VARIABLE SYSTEMS 

Entanglement in continuous variable systems has been extensively studied, see for instance p3. 12^ 12^ . We sum- 
marize here the results we will need in the remainder of the article. 

Consider two oscillators whose phase space variables {(Ir,pr) and {qL,PL) obey the canonical commutation relations. 
It is convenient to group the phase space variables as ^ = {qR,pR,qL,PL)- We can write the canonical commutation 
relations as 

where f2 is the symplectic matrix 






1 




























1 








-1 






(12) 
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For any quantum state p of the two oscillators, we can compute the first and second moments of its phase space 
variables 

< > = Trp^i (13) 
cri^ = ^Trp{^i- < >,^„r- < 6n >} (14) 

where {,} is the anticommutator. The covariance matrix of the oscillators is a real, symmetric, positive matrix, 
satisfying the constraint (which follows from positivity of the Hilbert Schmidt norm) 

a + i^Cl>0. (15) 

In general the first and second moments are only a partial characterisation of the quantum state p. But in the 
particular case where p is Gaussian, they completely characterise the state. 

The correlation matrix a allows one to study the entanglement of p. Denote by 



A 



PT 



/ 1 
10 
10 

Vo -1 



the matrix which realises the partial transpose. A necessary condition for entanglement of the two oscillators is 

ApTO-ApT + i^n < . (16) 

This condition is also sufficient if the oscillators are in a Gaussian state. It is convenient to rewrite this entanglement 
conditions as follows. Express the covariance matrix Hll|l as a bloc matrix 



then eq. (|16|l is equivalent (when A = B) to 



W = i(det[A] - det[C]) - det[a] - ^ < . (18) 

Below we will use the logarithmic negativity En as quantitative measure of entanglement. It is defined as 

i;Ar = max[0,-ln277"] (19) 

where rj^ is the smallest symplectic eigenvalues of the matrix KptO KpT '■ 

1 



S(t7) ± (I](a)2 - Adet{a)y/^ (20) 



with S(<t) — det{A) +det{Ji) — 2det{C) . The logarithmic negativity is an entanglement monotone. It is an upper bound 
on the distillable entanglement and a lower bound on the entanglement of formation. It measures the entanglement in 
units of entanglement bits (ebits), where one ebit is the entanglement present in a singlet state. Positivity of Ejsi > 
is a necessary condition for entanglement. In the particular case of Gaussian states it is both a necessary and sufficient 
condition. 

IV. CORRELATIONS BETWEEN TWO UNIFORMLY ACCELERATED OSCILLATORS 

To compute the elements of the covariance matrix first we have to quantize the quantum field in Rindler coordinates 
and evaluate several Minkowskian vacuum expectation values. We refer to [3, Q for detailed discussions of how to 
carry out these calculations, and summarize here very briefly the main points. 

To lighten the formulas we shall chose our length unit so that : 



a = 1 



(21) 
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The quantum field being massless, it decomposes into left (<— ) and right (— >) modes: 

$(c/,y) = + $ (y) . (22) 

These modes themselves split into modes defined on the left (L) and right (R) Rindler quadrants. For example : 

$(f/) = + . (23) 

We can decompose these operators in terms of Rindler modes, and Rindler creation and destruction operators 

r°° 1 

^^{u) = (e-^^-bU^) + e^'-bijX)) ^==dX , (24) 
Jo V47r A 

f°° 1 
^^{u) ^ {e^^-b^,{X) + e'^^-bl{X))^==dX . (25) 
^ Jo vAttX 

We emphasize the sign change in the arguments of the exponentials when we pass from the left to the right Rindler 
quadrant. It is the reflect of the opposite u-timc orientation in these quadrants. 

We can also decompose the field in terms of Minkowski modes and Minkowski creation and destruction operators 
^t/(v)H' bu{v)iwy. 

$([/)=/ (e-'"'^6t/(u;) + e+'"'%(u;)) ^=d«; . (26) 

The link between these two decompositions is provided by standard Bogoljubov transformations. Using the Bogoljubov 
transformations one shows that the Minkowski vacuum is perceived by the uniformly accelerated observer as being 
populated by a thermal bath of Rindler quanta at temperature Tjj = ^ . 

The Bogoljubov transformation also allows us to evaluate expectation values such as: 

^({*k(A),$«(A')}) = ({J^(A),J^(A')}) = ^cothA^<5(A + A') , (27) 

i({d«(A),$,(A')}> = ({J«W,!.(A')}> = ^-^<5(A + A') . (28) 

These expressions are then used to evaluate the correlation matrix between position and momentum variables expressed 
as in eqs. (UHl El) : 

, / ^ . ^^ 4 7 f+°° X^ COth AtT „ 

qqRB. = {qR{TR)qR{TR)) = / — ^— — ^— — — — rr dX (29) 

J-oo (A — i2+)(A — i2_)(A — i2+)(A — i2_) 4A7r 

(30) 

where we have introduced the notations: 
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2 



Amuj 



n± = uj{i-f± Vi-7^) , n± = cj{-i-f ± . (31) 



As shown in the appendix, the integral over A can be carried out, giving a closed form for qqRR. 
Similarly the other correlators can be expressed as 

qqR,{T) ^ {qR{rR)q,{r,)) ^ -—- - Cl^^X - 4n .inhXn '^^ ' ^^'^ 

1 ,r / N / 9 f^°° 1 COth Avr „ , , 

-^i{Mrn),PnirR)})=UWj_ (, _ ^,)(, _ ^_)(, _ ^,)(, _ 4. ^ ' ' ^''^ 

,Pn.iT) . -{{qUrMrR)}) (, . ^,)2(, . ^_). 4 . sinh a/ ^ ' ^''^ 

In the last two expressions there is a pole at A = 0. We resolve the ambiguity in the resulting integrals by integrating 
in the sense of a principal part, thereby obtaining for the first integral and a closed form given in the appendix for 
qPRL- 
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FIG. 2: Maximum entanglement between the oscillators, as a function of the dimensionless parameters 7 and lu. The outer line 
is the frontier of the region where the oscillators get entangled. The contour lines correspond to increases in the logarithmic 
negativity En in steps of 0.1. The maximum entanglement occurs for 7 ~ 0.703, ui ~ 0.0845 when En = 1.406. 



The computation of the momentum correlators are more delicate because they diverge due to the of a double pole 
at A = 0. We therefore introduce the infinite constant 

K^2^ f A-dA , (35) 



and so obtain: 

,,,,,, /■+00 I coth Att 

K+ppnn = {pn{rn)Pn{rn))-iimu; J^^ (a ^ (^^)(A - 1^_)(A - f^+)(A - »_) 4Xn ' ^^'^ 



K+ppnAT) ^ -i{pnirn),Pdr.)})^^,m.^l^ (A - r^,)^(A - ^^,)M A . sinh a/ ^ ' ^''^ 

where ppbr. and pprl (T) are finite quantities whose closed forms are given in the appendix. 
We may therefore write the correlation matrix as 

/ qqRR qqRL{T) qpRL{T) 

^ K + ppRR qpRL (T) K + PPRL {T) 

^ qqRL{T) qpRL{T) qqRR 
\ qpRL {T) K + PPRL [T) K + ppRR 

In this expression, except for K which is infinite, all the other terms are finite functions, depending only of the 
dimensionless parameters cj, 7, T, which can be interpreted as follows. The parameter oj is the ratio between the 
transition frequency of the oscillator and the acceleration a: u = ^ = "^here Tjj is the Unruh temperature. 

When cj > 27r the probability that the oscillator is excited will be exponentially small. We therefore expect any 
entanglement between the oscillators to disapear for large values of uj (since entanglement requires superpositions 
between several states). The parameter 7 = is the ratio between the line width F = (the inverse lifetime) of 
the first excited state of the oscillator and its transition frequency oj. When 7 > 1 the oscillator is strongly coupled to 
the field, whereas when 7 < 1 the oscillator is weakly coupled to the field. This translates in the Heisenberg equations 
of motion into the difference between the free solution being over damped or oscillating as it decays. In what follows 
we shall only consider the regime 7 < 1. Finally T = a{TR — tl) is the difference of the Rindler times along the 
two trajectories, in units of the inverse acceleration. It measures the Lorentz invariant distance of points on the two 
trajectories: As^ = 4 cosha T/2. 

We have checked numerically that the correlation matrix obeys the positivity constraint eq. H15|l . Indeed when 
inserting eq. H38|l into eq. (|15|l we find that the resulting matrix has one infinite positive eigenvalue, and three finite 
eigenvalues which we found to be positive using the procedure outlined in the appendix. 

Similarly we can consider the condition of positive partial transpose eq. (|16|l . Once more, inserting eq. (|38|l into 
eq. (|16() we find that the resulting matrix has one infinite positive eigenvalue, and three finite eigenvalues. When 
one of these eigenvalues becomes negative the state is entangled. We have also computed the logarithmic negativity 
En which quantifies the degree of entanglement present in the system. We find that logarithmic negativity is always 
finite and independent of K (see ea. l|Sn|) '). 
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FIG. 3: Evolution of the logarithmic negativity En between the two oscillators as a function of T — th ~ tl for the values 
7 = 0.703, uj — 0.0845 for which it reaches its maximal value of 1.406. In these units the period of the oscillator is 2Tr/u) ~ 74, 
and the "lifetime" of the first excited state of the oscillator is F"^ — {'yuj)~^ — 17. (We put quotes because the concept of 
lifetime is not well defined when it is smaller than the oscillation period, although it serves as a useful guiding concept for 
understanding what are the different time scales involved) . The oscillator is thus entangled for less than one oscillation period, 
and during this time it exchanges several quanta with the Unruh heat bath. Note also that in this case the Boltzman factor 
g-2Ti^ ~ 0.6 is of order 1 which implies that several of the oscillator levels have significant probability of being populated -this 
is necessary since En > 1 requires a system with dimension greater than 2-. The shapes of the curves for other values of 7 
and u) are similar, except that the maximum entanglements reached are smaller, and that the times at which these maxima are 
reached are different. 

Thus, even though the fluctuations of the oscillators coupled to the field are infinite, since the momentum correlators 
are infinite, the model is well defined. In particular the quantity we are interested in, the entanglement between the 
two oscillators, is always finite. 

We have computed the entanglement between the two oscillators as a function of T for different values of lu and 
7 < 1. We find that there are only specific pairs of values (7,0;) for which the detectors become entangled. In Fig. [21 
we have plotted the pairs of values of (7, w) for which the detector gets entangled. Note that our numerical analysis 
indicates that the region where entanglement occurs does not touch the axes 7 = and uj ~ 0. This is interesting 
since these axes correspond to the domain of validity of perturbation theory. Indeed the perturbative limit should 
arise when e — > 0, m fixed which corresponds to ^yuj — e^/m —^ 0. Thus the entanglement between two uniformly 
accelerated oscillators in opposite Rindler quadrants is a non perturbative phenomena. 

In Fig. 13 we have also plotted the degree of entanglement between the two oscillators as a function of (7,0;). We 
see that the entanglement is maximum for 7 ~ 0.703, lu ~ 0.0845 whereupon the logarithmic negativity reaches the 
value En = 1.406, see Fig. Ofor further discussion of this case. 

We have also computed how the entanglement between the two oscillators evolves as a function of T. For all values 
of the parameters, we find that the entanglement only appears when T > 0. That is the entanglement only gets 
established after a configuration of closest approach (T = tj^ — tl = 0) has been realized. The entanglement then 
increases, reaches a maximum, decreases and goes zero at a finite value of T. We can understand this as follows. 
The oscillators emit and absorb quanta that are packets localised in frequency and time (for instance if 7 << 1 then 
AA ~ £^/2m and At ~ (e^/2m)~^). These quanta are only correlated around configurations of closest approach 
T ~ 0. It takes some time to establish correlations between the two detectors, which is why the entanglement only 
appears for T > 0. At late times the quanta exchanged between the two oscillators are no longer correlated. The 
entanglement is gradually erased and finally disappears. As illustration in Fig. |3|we have plotted how the entanglement 
between the two detectors evolves as a function of T for the value of 7 and lu for which entanglement is maximal 

Note that because of boost invariance there is in fact an infinite set of pairs of locations for which the oscillators 
are entangled. Indeed if we look at one oscillator at a specific value of its proper time Ti(jj) = r*, it is only entangled 
with the other for a finite interval of proper time starting slightly after Tfl(L) = t* . This is depicted schematically in 

Fig. HI 

In summary we have studied the entanglement between two uniformly accelerated oscillators in 1+1 dimensional 
Minkowsky space time coupled to a massless scalar field. This model is exactly soluble. It allows us to study the 
case where both detectors are in thermal equilibrium with the field. It also allows us to study the case where the 
detectors are strongly coupled to the field. We find that there are some choices of parameters and of positions along 
the trajectories for which the two detectors get entangled. The maximum entanglement we find is slightly larger than 
1.4 entanglement bits. 
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V. APPENDIX: DETAILS OF CALCULATIONS 

We group here some of the calculations that are behind our main results. 

A. Explicit expressions for the correlation matrix elements 



Here we give explicit expressions for the correlators, which are expressed in terms of integrals in the main text. 
Specifically, upon carrying out the integrals in eqs. H29I32I34I36I37|I . one finds 



{qR{TR)qR{TB)) 



sinh ^27r-\/l — j^u 



27a; 



E 



n 



(cosh {l-K^X-i^ij^ - cos(27r7u;)) ^ „=i W' + ^'^ f - ^n^l^^ 



(39) 



(1 — 7"^) ' oj 
(-l)"ne-"^ 





(1 - 7^ + ?:7^1 - 72)w(r + ITT coth 


■nuj{^Jl — 7^ + 47) 


)-7) 


sinh 







27W 



(-l)"ne 



{m? + 27wn + w^)^ ^ 



(40) 



1 



7 0(T), 



2({ _ ^2)3/2 



3? 



JT^/T^u^ ^ ^1 _^2^(7r coth Tiuj{y/l -72 + 47) - iT)^ 



sinh |^7r(^l — 7^ + i^)uj 



27^3 



("^ -27^^" + ^^)^ 



(41) 



'<Pr{tr)pr{tr)) = A' + 



(^(1 - 272) sinh {211^1 - 72^) - 27^1 - 72 sin(27r7w)) 



2^/1 - 72 (cosh {21:^/1 - -I'^u^j - cos(27r7Cj)) 



27c^5g 



n=l n ((n2 + Ci)2^2 _ ^^2^2^2^ 



(42) 



^(WXrH),PL(r^)})=j^+ ^/^^^' 3,, 
2 (l-72)''/2 



3? 



87-2473 + 2(1-72)3/2 



sinh 



7r(-\/l — 72 + 17)0^ 



(1 — 72 — i 7-\/l — 72)ijJ ( T + iTT coth 7rti)(-\/l — 72 + 17) 



+ 



70;^ 



'[T] 2£ 



(-l)"e 



sinh 

47 — Tlj 



n (n2 — 2^ujn + lj2) 



7r(-\/l — 72 + i7)tjj 



n-T] 2^ 



(-l)"e 



n^nT 



(43) 



47 — Tw 



^^j^ 71 (n2 + 27wn + cj2) 



10 



B. Positivity and Entanglement 

Positivity of the Hilbert space inner product implies the positivity of the matrix a + iSl, see eq. H15|) . We checked 
that this is indeed the case as one of the eigenvalue of this matrix is infinite = while the three other are given by 
the eigenvalues of the matrix acting on the orthogonal space to the eigenvector of this infinite eigenvalue: 



qqRR^ - qqRL 5 ~ ^prl 

- 5 - qPRL PPRR ~ PPRL 














qqRL + qqRR 



(44) 



The three eigenvalues of this matrix are easy to compute. Using the above expressions for the correlators, we have 
checked numerically that they are positive, as expected. 

The criterium to put into evidence entanglement for a Gaussian system consists to show the occurence of negative 
eigenvalue in the partially transposed of the previous correlation matrix, i.e. the negativity of the matrix 



/ qqRR 
qpRR - i 
qqRL 
\ -qpRL 



qpRR + § 
K + ppRR 

qpRL 
-K - ppRL 



qqRL 
qpRL 
qqRR 
-qpRR - 



-qpRL \ 
-K - ppRL 
-qpRR + I 
K + ppRR J 



(45) 



Here again we find that one of the eigenvalues of this matrix is infinite = K. The computation of the three other 
eigenvalues is less obvious than in the previous case. First we may perform a symplectic transformation, using the 
matrix 



1 





1 


V2 







1 





1 





1 









1 

V2 






\ 

1 



71 / 



(46) 



to obtain the expression 

/ qqRR — qqRL 



2K 



-qpRL 




V 





— qpRL 



- 

2 " 

ppRL + ppRR qpRL 
qpRL qqRL + qqRR | 

-| ppRR-ppRL J 



(47) 



from which it is easy to isolate the eigenspace attached to the infinite eigenvalue and its orthogonal subspace. But to 
obtain the remaining three eigenvalues, i.e. the eigenvalues of the reduced matrix 



qqRR — qqRL ^qpRL 

V = I qqRL + qqRR f 

—qpRL — ^ ppRR — ppRL 



(48) 



we have to use (in principle) the general Cardan formula. We have performed such an analysis numerically. We have 
also used the criteria eq. ((T^ which in the present case reduces to 



lim R-'^W 

K^oc 



-2det[V] 



Finally the smallest symplectic eigenvalue of the partial transpose of a can be expressed as: 



V- 




dct[V] 



{qqRR - qqRL) 



(49) 



(50) 



It is independent of K. The expressions (|49|) and H50|l were used to compute numerically the results discussed in the 
main text. 
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